First Detection of Benign Rabbit Caliciviruses in Chile.

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Living with voracious roommates: Factors that explain isotopic niche variation in a mixed colony of insectivorous bats.

Theory predicts that in resource-limited environments, coexisting species may overlap their niche dimensions but must differ in at least one to avoid competitive exclusion. Specifically, it has been suggested that the coexistence of competing species within a guild, could be sustained with mechanisms of resource partitioning, such as segregation along a trophic dimension. Among the most gregarious mammals are bats, which present diversification in their diet based on habitat choice and body size. Despite differences that could explain specialization in prey selection, there are insufficient studies that explore food overlap in mixed bat colonies and the factors that determine the selection of prey, both at intra- and inter-specific levels. To fill this gap, we analyzed the isotope signal (δ13C and δ15N) in feces collected in a mixed colony of Tadarida brasiliensis and Myotis chiloensis. To understand how several factors could influence these isotopic signals, intrinsic explanatory variables were analyzed, including body mass, body length, age, and sex. Also, extrinsic variables were analyzed, including monthly temporality and moonlight intensity. Our findings support age-dependent specialization in M. chiloensis, with a significant role of moonlight intensity and sex on δ15N. In T. brasiliensis, we identified a significant effect of size, sex, and ear length on δ15N. Our analysis indicates that both species of bats experience diverse degrees of overlap through austral summer months, affected by several factors that explain the variability in their fecal isotopic signals.

Exploring the influence of density-dependence and weather on the spatial and temporal variation in common vole (Microtus arvalis) abundance in Castilla y León, NW Spain.

BACKGROUND: The common vole has invaded the agroecosystems of northwestern Spain, where outbreaks cause important crop damage and management costs. Little is yet known about the factors causing or modulating vole fluctuations. Here, we used 11 years of vole abundance monitoring data in 40 sites to study density-dependence and weather influence on vole dynamics. Our objective was to identify the population dynamics structure and determine whether there is direct or delayed density-dependence. An evaluation of climatic variables followed, to determine whether they influenced vole population peaks. RESULTS: First- and second-order outbreak dynamics were detected at 7 and 33 study sites, respectively, together with second-order variability in periodicity (2-3 to 4-5-year cycles). Vole population growth was explained by previous year abundance (mainly numbers in summer and spring) at 21 of the sites (52.5%), by weather variables at 11 sites (27.5%; precipitation or temperature in six and five sites, respectively), and by a combination of previous abundance and weather variables in eight sites (20%). CONCLUSIONS: We detected variability in vole spatiotemporal abundance dynamics, which differs in cyclicity and period. We also found regional variation in the relative importance of previous abundances and weather as factors modulating vole fluctuations. Most vole populations were cyclical, with variable periodicity across the region. Our study is a first step towards the development of predictive modeling, by disclosing relevant factors that might trigger vole outbreaks. It improves decision-making processes within integrated management dealing with mitigation of the agricultural impacts caused by voles. © 2023 Society of Chemical Industry.

European Rabbit Invasion in a Semi-Arid Ecosystem of Chile: How Relevant Is Its Role in Food Webs?

The European rabbit (Oryctolagus cuniculus) is one of the main invasive species in Chile, where it became naturalized ca. 150 years ago. Their high reproductive capacity, lack of specialist predators, and great adaptability favored the settlement of rabbits in diverse mainland and island ecosystems of the country. Recently, rabbits have become central players in semi-arid ecosystems, such as those represented in Las Chinchillas National Reserve in north-central Chile. We undertook to analyze the place and role of rabbits in the food web of that Reserve, based on a bibliographic review and long-term annual data gathered from 1987 until 2022 (36 years). Results showed that the network comprised 77 species, where 69% were primary producers (plants), 18% were mid-level consumers (herbivores), and 13% were top-level consumers (predators). The most connected species in the food web was the rabbit, which positively or negatively affected the species interacting with it. Predators such as Galictis cuja, Geranoaetus polyosoma, Leopardus colocolo, and Puma concolor, and the scavenger Vultur gryphus, could be negatively affected by an eventual decrease (natural or human-caused) in the rabbit population of the Reserve. To the contrary, primary producers such as Oxalis perdicaria, Plantago hispidula, Schizanthus parvulus, Senna cumminggi, and Tropaeolum azureum could be positively affected by an increase in their biomass in response to a decrease in rabbits, favoring native rodents. We consider that analyzing the rabbit-centered food web and its impacts on native interacting species allows a better understanding of the relevance of invasive species in the local community, providing conceptual tools for rabbit management.

Covariation of taxonomic and functional facets of ß-diversity in Chilean freshwater fish assemblages: Implications for current and future processes of biotic homogenization.

The biodiversity of assemblages that experience the introduction and extinction of species may lead to responses in two important facets: The taxonomic and functional diversity. The way in which these facets are associated may reveal important implications and consequences for the conservation of those assemblages. Considering the critical situation of freshwater fishes in continental Chile (30° - 56° S), we analyzed how the taxonomic (TDß) and functional (FDß) facets of ß-diversity, and their components of turnover and nestedness, are associated. We evaluated changes in ß-diversity (ΔTDß and ΔFDß), turnover (ΔTDtur and ΔFDtur), and nestedness (ΔTDnes and ΔFDnes) in 20 fish assemblages from their historical (pre-European) to current composition. We also simulated future trends of these changes, assuming that native species with conservation issues would become extinct. Our results show that the fish assemblages studied are in a process of loss of ß-diversity, both in taxonomic and functional facets (ΔTDß = -3.9%; ΔFDß = -30.4%); also, that these facets are positively correlated in the assemblages studied (r = 0.617; P < 0.05). Both components showed by loss in nestedness (ΔTDnes = -36.9%; ΔFDnes = -60.9%) but gain in turnover (ΔTDtur = 9.2%; ΔFDtur = 12.3%). The functional ß-diversity decreased more than the taxonomic (ΔFDß > ΔTDß), which was caused chiefly by six exotic species of Salmonidae, whose geographical spread was wider and that at the same time shared several morpho-functional traits. Our forecasts, assuming an intensification in the extinction of Endangered and Vulnerable native species, indicate that the process of homogenization will continue, though at a lower rate. Our study shows that the freshwater ichthyofauna of continental Chile is undergoing biotic homogenization, and that this process involves the facets of taxonomic and functional ß-diversity, which are show high correlation between historical and current compositions. Both facets show that process is influenced by nestedness, and while turnover contributes to differentiation (both taxonomic and functional), its importance is overshadowed by nestedness.

Partitioning ß-diversity reveals that invasions and extinctions promote the biotic homogenization of Chilean freshwater fish fauna.

AIM: Exotic species' introductions together with extinction of native species represent the main mechanisms driving biotic homogenization of freshwater fish assemblages around the world. While generally ichtyofaunistic realms transit towards biotic homogenization, for conservation purposes it is essential to understand what specific mechanisms are promoting it on particular areas or regions. Here, we report the occurrence of biotic homogenization in 29 Chilean watersheds, analyzing its ß-diversity (including turnover and nestedness) and predicting future trends. LOCATION: Continental Chile (18o-56o S). METHODS: We determined fish composition per basin for historical and current assemblages; extant native, exotic, and extinct species were recorded as 1 (presence) or 0 (absence) in two matrices basins × species. For each matrix, we calculated the turnover (ßsim), nestedness (ßnes), and ß-diversity (ßsor); then, we obtained Δßsim, Δßnes, and Δßsor, as the arithmetical difference between basin pairs over time. In addition, we search for explanatory variables correlating Δßsim, Δßnes, and Δßsor with geographical and land use variables. Finally, simulating events of species introduction (i.e., invasion) and extinction, we generated 15 hypothetical assemblages, looking to establish future trends towards biotic change in Chilean basins. RESULTS: Species turnover and ß-diversity significantly decreased from historical to current assemblages (Δßsim = -0.084; Δßsor = -0.061, respectively), while the species nestedness did not show significant changes (Δßnes = 0.08). Biotic changes have been driven mainly by the introduction of 28 exotic species, with a minor role of extinctions (one species) and translocations (0 species) of native species. Changes in ß-diversity were negatively correlated with area, elevation, and geographical distance between basins but not with land-use nor human population. Finally, the analysis of 15 future assemblages predicts a significant decrease of ß-diversity and turnover, and an increase for species nestedness, this time promoted by an increase in the extinction of native species. MAIN CONCLUSION: Chilean basins show a significant decrease of the distributional ß-diversity and species turnover of the freshwater fish fauna, evidencing a trend towards biotic homogenization. This trend is shared with other Neotropical basins; however, specific mechanisms driving it show different magnitude. Changes in the ß-diversity components do not show correlation with variables associated to land use, thus suggesting that casual introductions of freshwater fishes in Chile follow an opportunistic mode related to commercial use. According to future scenarios simulated, biotic homogenization should increase further, mainly as consequence of increased native extinctions.

Latitudinal patterns in the diet of Andean condor (Vultur gryphus) in Chile: Contrasting environments influencing feeding behavior.

Human-dominated environments alter the availability and quality of resources for many species, especially for scavengers that have large home ranges and plastic foraging behaviors that enable them to exploit novel resources. Along the western slope of the Andes, the modification of natural landscapes have resulted in significant declines in native prey, the introduction of non-native species, and an increase in the availability of anthropogenic resources. These factors have likely influenced the resources available to Andean condors (Vultur gryphus), however, data are lacking as to how condor's diet vary along their large latitudinal range. We evaluated differences in Andean condor diet along a ~2500 km latitudinal gradient in Chile from the heavily modified Central zone (32-34°S) to the more pristine Austral zone (44-56°S). We assessed diet composition through the identification of prey remains in condor pellets, and carbon and nitrogen isotope analysis of condor feathers and their primary prey identified from pellet analysis. Our results identified medium- and large-bodied domesticated mammals (ungulates) and introduced exotic species (lagomorphs) as common prey across the study area. Condors from the Central zone had the largest isotopic niche width, probably related to consumption of anthropogenic resources with distinctly high carbon isotope values indicative of C4-based foods likely acquired from landfills or corn-fed livestock. Isotopic niches for condors from the Southern and Austral zones almost completely overlapped. Andean condor diet is strongly influenced by local conditions determining differential access to prey sources. The high dependence of Andean condors on livestock across a large geographical area, and landfills in more (sub)urban areas, may help stabilize their populations via anthropogenic resources subsidies. Long-term dependence on such resources, however, may have health costs including contaminant exposure and greater mortality risk. These data will help identify potential threats related to resource availability and use, and better inform management and conservation decisions.

Extreme climatic events change the dynamics and invasibility of semi-arid annual plant communities.

Extreme climatic events represent disturbances that change the availability of resources. We studied their effects on annual plant assemblages in a semi-arid ecosystem in north-central Chile. We analysed 130 years of precipitation data using generalised extreme-value distribution to determine extreme events, and multivariate techniques to analyse 20 years of plant cover data of 34 native and 11 exotic species. Extreme drought resets the dynamics of the system and renders it susceptible to invasion. On the other hand, by favouring native annuals, moderately wet events change species composition and allow the community to be resilient to extreme drought. The probability of extreme drought has doubled over the last 50 years. Therefore, investigations on the interaction of climate change and biological invasions are relevant to determine the potential for future effects on the dynamics of semi-arid annual plant communities.

Low functional richness and redundancy of a predator assemblage in native forest fragments of Chiloe island, Chile.

1. Changes in land use and habitat fragmentation are major drivers of global change, and studying their effects on biodiversity constitutes a major research programme. However, biodiversity is a multifaceted concept, with a functional component linking species richness to ecosystem function. Currently, the interaction between functional and taxonomic components of biodiversity under realistic scenarios of habitat degradation is poorly understood. 2. The expected functional richness (FR)-species richness relationship (FRSR) is positive, and attenuated for functional redundancy in species-rich assemblages. Further, environmental filters are expected to flatten that association by sorting species with similar traits. Thus, analysing FRSR can inform about the response of biodiversity to environmental gradients and habitat fragmentation, and its expected functional consequences. 3. Top predators affect ecosystem functioning through prey consumption and are particularly vulnerable to changes in land use and habitat fragmentation, being good indicators of ecosystem health and suitable models for assessing the effects of habitat fragmentation on their FR. 4. Thus, this study analyses the functional redundancy of a vertebrate predator assemblage at temperate forest fragments in a rural landscape of Chiloe island (Chile), testing the existence of environmental filters by contrasting an empirically derived FRSR against those predicted from null models, and testing the association between biodiversity components and the structure of forest fragments. 5. Overall, contrasts against null models indicate that regional factors determine low levels of FR and redundancy for the vertebrate predator assemblage studied, while recorded linear FRSR indicates proportional responses of the two biodiversity components to the structure of forest fragments. Further, most species were positively associated with either fragment size or shape complexity, which are highly correlated. This, and the absence of ecological filters at the single-fragment scale, rendered taxonomically and functionally richer predator assemblages at large complex-shaped fragments. 6. These results predict strong effects of deforestation on both components of biodiversity, potentially affecting the functioning of remnants of native temperate forest ecosystems. Thus, the present study assesses general responses of functional and taxonomic components of biodiversity to a specific human-driven process.

Effect of prolonged exposition to pulp mill effluents on the invasive aquatic plant Egeria densa and other primary producers: a mesocosm approach.

The recent disappearance of the aquatic plant Egeria densa, a Brazilian native invasive species, from a wetland in southern Chile prompted several efforts to unveil the origin of this phenomenon. Because these changes occurred by the time a newly built pulp mill started its operations in the area, a reasonable doubt for a cause-effect relationship is plausible. We implemented a mesocosm approach to directly evaluate the effect of treated pulp mill effluents (PMEs) on several growth-related parameters of E. densa as well as other primary producers. We hypothesize that effluent, at a dilution similar to that detected in the zone of the wetland where the negative environmental impacts were evident, has a significant negative effect on this aquatic plant as well as on other primary producers inside a mesocosm system. After a prolonged (months) exposure to both 0% PME with pure river water and a 4 to 5% (v/v) dilution of treated PME, no effect on E. densa was measured. Furthermore, plants exposed to effluent exhibited a significantly greater general growth rate. Coincidently, chlorophyll a concentration in the water column and periphyton biomass also changed over time, but without any pattern attributable to the effluent. Values of the autotrophic index obtained from the periphyton growth pattern did not suggest enrichment of the system with organic matter. Our results only refer to the direct effect of mill effluents on several biotic responses, but they represent an important advance toward generation of the scientific knowledge necessary to understand how the ecosystem functions while receiving this and other unquantified sources of water.

Effects of functional constraints and opportunism on the functional structure of a vertebrate predator assemblage.

1. Within mainstream ecological literature, functional structure has been viewed as resulting from the interplay of species interactions, resource levels and environmental variability. Classical models state that interspecific competition generates species segregation and guild formation in stable saturated environments, whereas opportunism causes species aggregation on abundant resources in variable unsaturated situations. 2. Nevertheless, intrinsic functional constraints may result in species-specific differences in resource-use capabilities. This could force some degree of functional structure without assuming other putative causes. However, the influence of such constraints has rarely been tested, and their relative contribution to observed patterns has not been quantified. 3. We used a multiple null-model approach to quantify the magnitude and direction (non-random aggregation or divergence) of the functional structure of a vertebrate predator assemblage exposed to variable prey abundance over an 18-year period. Observed trends were contrasted with predictions from null-models designed in an orthogonal fashion to account independently for the effects of functional constraints and opportunism. Subsequently, the unexplained variation was regressed against environmental variables to search for evidence of interspecific competition. 4. Overall, null-models accounting for functional constraints showed the best fit to the observed data, and suggested an effect of this factor in modulating predator opportunistic responses. However, regression models on residual variation indicated that such an effect was dependent on both total and relative abundance of principal (small mammals) and alternative (arthropods, birds, reptiles) prey categories. 5. In addition, no clear evidence for interspecific competition was found, but differential delays in predator functional responses could explain some of the unaccounted variation. Thus, we call for caution when interpreting empirical data in the context of classical models assuming synchronous responses of consumers to resource levels.

On the relationship between productivity and food chain length at different ecological levels.

The effects of energy on food web structure have been debated for at least 80 years. Nevertheless, the empirical evidence is meager, especially from terrestrial ecosystems. We analyzed long-term temporal variation in food chain length in a semiarid continental ecosystem, where productivity shows large interannual variations. Incidence of nonherbivorous prey in predator diet was used as a proxy of trophic position, allowing us to analyze the effect of productivity on food chain length within the assemblage of top predators (which comprises the most abundant and persistent top predators in the system) and to compare observed patterns at the species and assemblage levels. At the species level, the relationship between trophic position and productivity took different forms, varying in magnitude and shape. This pattern contrasts with the consistent increase in food chain length, with productivity observed at the assemblage level. Our results indicate that productivity can be a main determinant of food chain length, but not necessarily because of energy limitation. Further, the increase in food chain length with available energy probably represents an aggregate attribute, driven to a large extent by predators with higher consumption rates, rather than being the result of compensatory responses among predators.

Interplay between metabolic rate and diet quality in the South American fox, Pseudalopex culpaeus.

We studied the metabolic costs associated with the ingestion of peppertree fruits (Schinus molle) in the culpeo fox, Pseudalopex culpaeus, the second largest canid in South America. Throughout its range of distribution, this fox feeds on rodents and other small vertebrates, and also on peppertree fruits, which represent 98% of total fruits consumed in semiarid Chile. Peppertree contains a high diversity of phytochemicals. Foxes feeding on diets containing rats and peppertree fruits (mixed diets) exhibited a 98.9% increase in basal rate of metabolism when compared to rat-acclimated foxes. Thus, acute ingestion of chemically defended fruits has an energetic cost for the fox, reflected in higher values of basal metabolism. Increased metabolic rates may be associated with increased protein synthesis for detoxification and for tissue repair, including the production of biotransformation enzymes.

Population dynamics of small mammals in semi-arid regions: a comparative study of demographic variability in two rodent species.

The seasonally determined demographic structure of two semi-arid rodents, both agricultural pest species (the leaf-eared mouse (Phyllotis darwini) in Chile and the multimammate mouse (Mastomys natalensis) in Tanzania), is analysed using capture-mark-recapture (CMR) statistical models and measures for elasticity (the relative change in the growth rate due to a relative unit change in the parameter of concern) derived from projection linear matrix models. We demonstrate that reproduction and survival during the breeding season contribute approximately equally to population growth in the leaf-eared mouse, whereas the multimammate mouse is characterized by a more clearly defined seasonal structure into breeding and non-breeding seasons and that reproduction contributes far more than survival during the breeding season. On this basis, we discuss evolutionary and applied (pest control) issues. Regarding the evolution of life histories (leading to a maximization of the overall net annual growth rate), we suggest that for the leaf-eared mouse, features favouring survival throughout the year will provide selective value, but that during the main breeding season, features favouring reproduction and survival are about equally favourable. For the multimammate mouse, features favouring survival are particularly important outside the breeding season, whereas during the breeding season features favouring reproduction are more important. Regarding pest control (aiming at reducing the overall net annual growth rate), we suggest that (ignoring economic considerations) affecting survival outside the main breeding season is particularly effective for the leaf-eared mouse, a feature that is even more the case for the multimammate mouse. In sum, we demonstrate through this comparative study that much is to be learnt from studying the dynamics of fluctuating small rodents-a focal issue within much of population ecology.

Population dynamics of a South American rodent: seasonal structure interacting with climate, density dependence and predator effects.

Understanding the role of interactions between intrinsic feedback loops and external climatic forces is one of the central challenges within the field of population ecology. For rodent dynamics, the seasonal structure of the environment necessitates changes between two stages: reproductive and non-reproductive. Nevertheless, the interactions between seasonality, climate, density dependence and predators have been generally ignored. We demonstrate that direct climate effects, the nonlinear effect of predators and the nonlinear first-order feedback embedded in a seasonal structure are key elements underlying the large and irregular fluctuations in population numbers exhibited by a small rodent in a semi-arid region of central Chile. We found that factors influencing population growth rates clearly differ between breeding and non-breeding seasons. In addition, we detected nonlinear density dependencies as well as nonlinear and differential effects of generalist and specialist predators. Recent climatic changes may account for dramatic perturbations of the rodent's population dynamics. Changes in the predator guild induced by climate are likely to result, through the food web, in a large impact on small rodent demography and population dynamics. Assuming such interactions to be typical of ecological systems, we conclude that appropriate predictions of the ecological consequences of climate change will depend on having an in-depth understanding of the community-weather system.

El Niño-Southern Oscillation-Driven Rainfall Variability and Delayed Density Dependence Cause Rodent Outbreaks in Western South America: Linking Demography and Population Dynamics.

It is well known that some rodent populations display dramatic density fluctuations in semiarid regions of western South America after the unusual rainfall levels associated with El Niño-southern oscillation (ENSO) disturbances. These correlated phenomena have led some ecologists to believe that rodent outbreaks are determined solely by density-independent factors (e.g., rainfall regime). However, demographic studies have detected strong delayed density-dependent effects in one of the most irruptive rodent species, the leaf-eared mouse Phyllotis darwini. We tested the effects of rainfall and delayed density-dependent factors by constructing a structured model based on demographic data estimated from a capture-mark-recapture study of this species in Chile. A model including both rainfall and delayed density-dependent effects predicts the observed population dynamics rather accurately over a 10-yr period. Interestingly, small changes in model parameters result in large changes in model dynamics, which strongly suggests that local variations in demographic features are important in explaining the asynchronous pattern in outbreak occurrences. These findings suggest that inextricably intertwined endogenous and exogenous forces cause rodent outbreaks in western South America. The former are characterized by delayed nonlinear feedbacks, whereas the latter are characterized by the positive effects of the El Niño phases and the negative effects of the La Niña phases of the ENSO disturbance.

Population extinction risks of three Neotropical small mammal species.

The population persistence and extinction probabilities of three small mammal species were analyzed by estimating growth and extinction properties obtained from 10 years of live-trapping data at two different habitat types in semiarid Chile. We used a stochastic formulation with an exponential growth model known as a Wiener-drift process, out of which growth and extinction quantities were estimated. The rodent Phyllotis darwini showed the lowest rates of growth, and the lowest infinitesimal variances, whereas the opposite trend was found for the rodent Akodonolivaceus. The marsupial Thylamys elegans showed intermediate values for growth rates and infinitesimal variances. The rodent P. darwini showed the lowest extinction risk in the study site. We also detected spatial differences between mesic and xeric habitats in the growth rates of P. darwini and T. elegans, and in the extinction risks of the three species studied. Although the population growth of these three species can be approximated by purely stochastic processes, the introduction of density dependence through autoregressive log-linear models reduced the extinction times of all species analyzed.

Extinction and colonization processes in subpopulations of five neotropical small mammal species.

We report the extinction and colonization rates of five sympatric small mammal species at a semiarid locality in north central Chile. We provide information based on 6 years of monitoring on how colonization and extinction rates change according to landscape features (slope aspect) and on their relationship to populations size, population variability, and body size. We found that: (1) for all species in the assemblage, extinction rates of subpopulations from equatorial-facing slopes were significantly lower than those in polar-facing slopes, (2) population size was the most important factor determining extinction rates, (3) colonization rates did not vary between slopes, and were affected by population size only in equatorial-facing slopes, and (4) most species had higher extinction than colonization rates. Persistence of the metapopulation system for all five small mammal species appears to be fueled by repeated colonization events.